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We give an example of a purely bosonic model - a rotor model on the 3D cubic lattice - whose low 
energy excitations behave like massless (7(1) gauge bosons and massless Dirac fermions. This model 
can be viewed as a "quantum ether": a medium that gives rise to both photons and electrons. It 
illustrates a general mechanism for the emergence of gauge bosons and fermions known as "string- 
net condensation." Other, more complex, string-net condensed models can have excitations that 
behave like gluons, quarks and other particles in the standard model. This suggests that photons, 
electrons and other elementary particles may have a unified origin: string-net condensation in our 
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I. INTRODUCTION 



Throughout history, people have attempted to under- 
stand the universe by dividing matter into smaller and 
smaller pieces. This approach has proven extremely fruit- 
ful: successively smaller distance scales have revealed 
successively simpler and more fundamental structures. 
Over the last century, the fundamental building blocks of 
nature have been reduced from atoms to electrons, pro- 
tons and neutrons, to most recently, the "elementary" 
particles that make up the 17(1) x SU(2) x SU(3) stan- 
dard model. Today, a great deal of research is devoted 
to finding even more fundamental building blocks - such 
as superstrings. 

This entire approach is based on the idea of reduction- 
ism - the idea that the fundamental nature of particles 
is revealed by dividing them into smaller pieces. But 
reductionism is not always useful or appropriate. For ex- 
ample, in condensed matter physics there are particles, 
such as phonons, that are collective excitations involv- 
ing many atoms. These particles are "emergent phenom- 
ena" that cannot be meaningfully divided into smaller 
pieces. Instead, we understand them by finding the mech- 
anism that is responsible for their emergence. In the case 
of phonons, for example, this mechanism is symmetry 
breaking. [1-4] 

This suggests alternate line of inquiry. Could the el- 
ementary particles in the standard model be analogous 
to phonons? That is, could they be collective modes of 
some "structure" that we mistake for empty space? 

Recent work suggests that they might be. [5, 6] This 
work has revealed the existence of new and exotic phases 
of matter whose collective excitations are gauge bosons 
and fermions. The microscopic degrees of freedom in 
these models are spins on a lattice - purely local, bosonic 
objects with local interactions. There is no trace of gauge 
boson or fermion degrees of freedom in the underlying 
lattice model. The gauge bosons and fermions are thus 
emergent phenomena - a result of the collective behavior 
of many spins. 
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What is the mechanism responsible for their emer- 
gence? In these exotic phases, the spins organize into 
a special pattern - a particular kind of entangled ground 
state, which we call a "string-net condensed" state. A 
string-net condensed state is a spin state where the spins 
organize into large string-like objects (or more generally 
networks of strings). The strings then form a quantum 
string liquid (see Fig. 1). This kind of ground state natu- 
rally gives rise to gauge bosons and fermions. The gauge 
bosons correspond to fluctuations in the strings - the col- 
lective motions of the strings that fill the space. [7-10] 
The fermions correspond to endpoints of the strings - 
that is, defects in the string liquid where a string ends in 
empty space. [5, 6, 11] 

What makes the string-net picture particularly com- 
pelling is that the gauge bosons and fermions naturally 
emerge together. They are just different aspects of the 
same underlying structure. Therefore, if we believe that 
the vacuum is such a string-net condensate then the pres- 
ence of gauge interactions and Fermi statistics in the 
standard model is no longer mysterious. String- net con- 
densation explains what gauge bosons and fermions are, 
why they exist, and why they appear together. [6] 

The general theory of string-net condensation was 
worked out in Rcf. [5]. One of the main results in that 
paper was a series of exactly soluble models realizing all 
possible string-net condensates. These models are quite 
general and can realize gauge bosons with any gauge 
group. However, they are also complicated when dis- 
cussed in full generality, and Ref. [5] did not provide an 
explicit example of the most physically relevant case - 
a model realizing gauge bosons and fermions in (3+1) 
dimensions. 

In this paper, we attempt to remedy this problem. We 
demonstrate the string-net picture of (3+l)D emerging 
gauge bosons and fermions with concrete lattice mod- 
els. We describe a rotor model on the cubic lattice that 
produces both i7(l) gauge bosons and fermions. The 
fermions can be gapped excitations (as in an insulator) 
or gapless (as in a Fermi liquid). They can also behave 
like massless Dirac fermions. In this case, the low energy 
physics of the rotor model is identical to massless quan- 
tum electrodynamics (QED). The rotor model can then 
be viewed as a "quantum ether": a medium that gives 
rise to both photons and electrons. In addition, the ro- 
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FIG. 1: A picture of quantum ether: the fluctuations of the 
string-nets give rise to gauge bosons (such as photons). The 
ends of the strings give rise to fermions (such as electrons) . 



tor model is closely related to U(l) lattice gauge theory 
coupled to a Higgs field. It demonstrates that a simple 
modification or "twist" can change the Higgs boson into 
a fermion. 

While this is not the first lattice bosonic model with 
emergent massless gauge bosons and massless Dirac 
fermions [7, 12-14], it has two characteristics which dis- 
tinguish it from previous examples. First, the mapping 
between the rotor model and QED is essentially exact, 
and does not require a large N limit or similar approx- 
imation. Second, the rotor model is a special case of 
a general construction, [5, 6], unlike the other models 
which were in some sense, discovered by accident. It 
therefore provides a deeper understanding of emergent 
fermions and gauge bosons. In addition to its relevance to 
high energy physics, this understanding may prove useful 
to condensed matter physics, particularly the search for 
phases of matter with exotic low energy behavior. 

The paper is organized as follows: we begin with a 
"warm-up" calculation in section II - a rotor model with 
emergent photons and bosonic charges. This model is 
closely related to U(l) lattice gauge theory coupled to a 
Higgs field. Then, in section III, we show that the ro- 
tor model can be modified in a natural way by adding 
a phase factor or "twist" to a term in the Hamiltonian. 
This modified or "twisted" rotor model has (massive) 
fermionic charges. In section IV, we describe further 
modifications which give rise to gapless Fermi liquids and 
massless Dirac fermions. 
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FIG. 2: A picture of the rotor model (1). The term Qi = 
( — I) 1 X^iegs of l^l ac ts on the six "legs" of I - that is, the 
six rotors adjacent to I, drawn above as filled dots. The term 
acts on the four rotors, labeled by 1, 2, 3, 4, 
along the boundary of the plaquette p. 




FIG. 3: A picture of a closed string state. The rotors along the 
thick closed curve have angular momentum L z — ±1, while 
all the other rotors (depicted as empty dots) have L z — 0. 



A. The model 

The model is a quantum rotor model with rotors on 
the links of a cubic lattice. Each rotor can be viewed as 
a particle moving on a circle. The position of the particle 
is given by the angle 9, and the angular momentum of 
the particle by L z = 19, where I is the moment of inertia. 
The Hamiltonian of the rotor model is given by (see Fig. 
2a) [8, 15] 



Bp — L 2 L 3 L 4 , 



legs of I 



II. A 3D ROTOR MODEL WITH EMERGENT 
PHOTONS AND BOSONIC CHARGES 



where I = [I x , I y ,I z ) labels the vertices, i labels the links 
and p labels the plaquettes of the cubic lattice. The "legs 
of J" are the six links that are attached to the vertex I, 
and 1, 2, 3, 4 label the four links of the plaquette p (see 
Fig. 2). L ± = e 10 are the raising/lowering operators of 
L z and = (-l)*«+'v+-r*. 



In this section, we present a warm-up example - a ro- 
tor model with emergent photons and bosonic charges. 
This model is closely related to U(l) lattice gauge the- 
ory coupled to a Higgs field. While a number of other 
similar models [8, 15-18] have been analyzed previously, 
this model has the advantage of being quasi-exactly sol- 
uble, and generalizing easily to the fermionic case. The 
string-net picture of the ground state is also particularly 
evident in this case. 



B. The string- net picture 

In this section, we show that the rotor model exhibits 
string- net condensation in the regime V ^> g ^> J > 0. 
In particular, we show that the rotors organize into ef- 
fective extended objects, namely "string-nets", and the 
ground state is a quantum liquid of these string-nets. The 
string-net condensed ground state has two types of exci- 
tations - string collective motions (which will correspond 
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to photons) and string endpoints (which will correspond 
to charges). 

The first step is to notice that the Qi operators com- 
mute with each other as well as the other terms in the 
Hamiltonian. Thus, we can label all the eigenstates of 
#rotor by their eigenvalues under Qf. {qi}. The quan- 
tum number qi - which we will call the "charge" on site 
I - takes values in the integers. 

Clearly the lowest energy charge configuration is the 
charge-0 configuration qi = 0. Other charge configura- 
tions cost an additional energy V. Hence, in the limit 
V ^> J, g, the low energy physics is completely contained 
in the charge-0 sector. 

We would like to enumerate the states in the charge-0 
sector. It is natural to describe these states in terms of 
extended objects - in particular, string-nets. The sim- 
plest state has Lf = for every rotor. In the string-net 
language, we think of this state as the vacuum. Another 
state can be obtained by alternately increasing or de- 
creasing L z by 1 along a loop (see Fig. 3). We think of 
this state as the vacuum with a single (closed) string. 
Other states can be constructed by repeating this pro- 
cess. One finds that the most general charge-0 state can 
have many strings and that the strings can overlap to 
form networks of strings or "string-nets." 

What is the action of the Hamiltonian (1) on these 
charge-0 states? If we think about the states using the 
string-net language, we see that the J-term penalizes 
strings for being long. It therefore corresponds to a string 
tension. On the other hand, the effect of the B p oper- 
ator is to either create a small loop of small string (if 
applied to the vacuum) or to deform the existing strings 
(if applied to a more complicated state). Thus the g- 
tcrm generates string "hopping" or string fluctuations. 
One can think of it as a string kinetic energy. 

There are two regimes to consider. When J ^> g, the 
string tension dominates and the ground state contains 
almost no strings. The ground state is a "normal" state. 
When g ^> J, the string kinetic energy dominates and 
the ground state is a superposition of many large closed 
strings [7, 8]. The ground state is thus "string-net con- 
densed." We expect a quantum phase transition between 
these two states at some J/g of order 1. 

We expect that the excitations above the ground state 
have very different properties in these two regimes. In 
the string-net condensed case, there are two types of ex- 
citations: low energy excitations with charge-0 and high 
energy excitations with nonzero charge. 

The low energy excitations can be constructed from 
linear combinations of string-net states. These excita- 
tions can be thought of as collective motions of closed 
strings. 

The high energy, charged, excitations can be con- 
structed from linear combinations of string-net config- 
urations with open strings. For example, if one takes 
the vacuum state L z = and alternately increases and 
decreases L z along an open path P, the resulting state 
has nonzero charge at the two ends of P. This is true 
quite generally: nonzero charge configurations are made 
up of open string-nets - with the charge located at the 
endpoints of the open strings. 



We can think of the string endpoints as point defects in 
the condensate. The higher energy excitations are com- 
posed out of these point defects. The defects carry a 
quantum number - charge - which is measured by the 
operator Qi. 

C. Equivalence with 17(1) gauge theory 

In this section, we examine the rotor model more quan- 
titatively. We show that the low energy physics of the 
rotor model (1) is described by compact U(l) gauge the- 
ory coupled to infinitely massive charges. The two types 
of excitations discussed in the previous section have a 
natural interpretation in terms of the gauge theory. The 
collective motions of the strings give rise to gapless exci- 
tations and behave like photons, while the point defects 
(the endpoints of the strings) are gapped and behave like 
charged particles. 

In the gauge theory language, the strength of the 
interaction between the photons (the gauge field) and 
the charges is characterized by a "fine structure con- 
stant" a. The "fine structure constant" also characterizes 
the strength of quantum fluctuations of the gauge field. 
When a>l, the quantum fluctuations are so large that 
the gauge theory is in the confining phase and there are 
no gauge bosons at low energies. The gauge bosons (such 
as photons) can exist at low energies only in the decon- 
fined phase when a is small. 

In the following, we will show that the "fine struc- 
ture constant" a of the gauge theory is of order \J J/g. 
The "speed of light" c is of order \J gJa? where a is the 
lattice constant. This result provides insight into the 
phase transition between the normal and string-net con- 
densed states: the normal state corresponds to the large 
a confining phase of the gauge theory, while the string- 
net condensed state corresponds to the small a decon- 
fined phase. The string-net condensation transition is the 
usual confinement-deconfincmcnt transition from gauge 
theory 

The simplest way to derive these results is to map (I) 
onto a lattice gauge theory Hamiltonian. Consider the 
following Hamiltonian for compact U(l) gauge field cou- 
pled to a charge- 1 Higgs field e 1 ^: 

i (U) 
- 2g cos(A IJ + Ajk + Akl + Ali) 

{IJKL) 

-- 2t^2cos((j) I -<Pj-A I j) (2) 
(U) 

Here, (IJ), (IJKL) is an alternative notation for links 
i, and plaquettes p in the cubic lattice. The operator Ejj 
is the integer valued electric field, canonically conjugate 
to the vector potential: [Au,Eu] = i. Similarly, nj 
is the occupation number operator canonically conjugate 
to the Higgs field: [</>/, nj] = i. 

In lattice gauge theory, one is interested in the proper- 
ties of this Hamiltonian within the gauge invariant sub- 
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space where Gauss' law holds: 

J2 E u = n J ( 3 ) 

(U) 

Our claim is that within this gauge invariant subspace, 
-Hgauge is mathematically equivalent to H I0t0 i- To see 
this, note that the electric field operators and the occu- 
pation number operators all commute with each other. 
So the lattice gauge theory model has a complete basis 
consisting of states |{e/j, nj}) with electric field e/j, and 
occupation number nj, 

Eij\{eij,m}) = e/j|{e7j,n 7 }) 
ni\{eu, nj}) = n 7 |{ejj,nj}) 

Similarly, the rotor model H rotm has a complete basis 
consisting of states |{/|}) with angular momentum If. 

On can map the rotor model Hilbert space onto the 
gauge invariant subspace of the lattice gauge theory 
by mapping the basis state onto the basis state 

|{e/j,n/}) given by 

eu = -If 

ni = E l l ( 4 ) 

legs of I 

(Here i is the link connecting sites J and J) . Clearly, this 
correspondence maps the operators nj and Eu onto the 
operators Qz and (— l) 1 ■ Lf. One can also show that 
this correspondence maps cos(Au + Ajk + Akl + Aei) 
onto l/2(B p + h.c.) and cos(0/ — <f)j — Ajj) onto cos(0j). 
Examining the two Hamiltonians, we conclude that the 
gauge theory in the special case t = can be mapped 
onto the rotor model. 

The physics of the rotor model (1) is therefore com- 
pletely equivalent to a U(l) gauge field coupled to 
bosonic charges. The fine structure constant a and speed 
of light c can be computed using the lattice gauge theory 
Hamiltonian (2)). [7, 8] As for the mass of the charges, 
the rotor model (1) corresponds to the special case where 
t = - the case where the charges are infinitely massive 
and have no dynamics. 

D. Adding a charge hopping term 

In the previous section, we discovered that the point 
defect excitations have an infinite mass. This is not sur- 
prising, since the operator Qj commutes with the Hamil- 
tonian (1), and hence charge configurations are com- 
pletely static. 

We would like to modify H xoioI so that the point defects 
(or "charges") do have dynamics. To do this, we need 
to add a term, i?ho P , to the Hamiltonian that doesn't 
commute with Qj. Such a term will make the charges 
hop from site to site. We expect that any i?h op will give 
rise to the small qualitative behavior (e.g. same particle 
statistics). We therefore consider the simplest hopping 
term: 

ffhop = -2tJ2co S (6i) = -tJ2(Lt + LT) 

i i 



It's easy to see that this term makes the charges hop 
from site to site just as in a nearest neighbor tight 
binding model. In fact, if we examine the mapping, 
(4), we see that Hh op corresponds exactly with the 
nearest neighbor hopping term in lattice gauge theory, 
cos(Afe0(r) — Ak{r)). Therefore, the low energy physics 
of the rotor Hamiltonian (H Iotm + -ffhop) is equivalent to 
that of a U(l) gauge field coupled to charges with finite 
mass. 



E. The statistics of the (7(1) charges 

In this section, we compute the statistics of the charges 
and show that they are bosons. The main purpose of 
this computation is to facilitate comparison with the 
"twisted" rotor model. (Indeed, the fact that the charges 
are bosons follows immediately from the exact mapping 
to lattice gauge theory, (4)). 

It is easiest to perform the computation in the case 
J = 0. In this case, B p commutes with the Hamiltonian 
(iJ rotor + -ffhop)- We can therefore divide the Hilbert 
space into different sectors {9 p } corresponding to differ- 
ent flux configurations B p = e t6p . A complete basis for 
each sector can be obtained by listing all the different 
charge configurations {qi}. The action of the Hamil- 
tonian (i/ r otor + -Shop) on a charge configuration {qi} is 
simple. The first part of the Hamiltonian, H rot0 r, doesn't 
affect the charge configuration at all, while the second 
part, -f/hop acts in two ways: it either creates two charges 
at neighboring sites, or it makes a charge hop from one 
site to another. 

Thus, within each sector, (i? ro tor + #hop) is simply a 
hopping Hamiltonian on the cubic lattice. The hopping 
operators L i make the charges hop from site I to site J, 
where I, J are the two endpoints of i. To compute the 
statistics of the charges, we use the statistical hopping 
operator algebra described in Ref. [11]. We note that the 
hopping operators satisfy the algebra 

L+LJL+ = L+LjLt (5) 

for any i,j,fc incident to some vertex I. According to 
Ref. [11], this is the bosonic hopping operator algebra. 
We conclude that the charged particles are bosons. 

III. A 3D ROTOR MODEL WITH EMERGENT 
PHOTONS AND FERMIONIC CHARGES 

In this section, we present another rotor model which 
is very closely related to the model from the previous 
section, (1). The Hamiltonian only differs from -ff ro tor by 
a "twist" - an additional phase factor in the B p term. 
Despite the apparent similarity, we will see that this 
"twisted" rotor model gives rise to emergent photons and 
fermionic charges. 

The relationship between the twisted and untwisted 
models is a special case of a general relationship between 
fermionic and bosonic string-net condensates. [5, 6] It 
is part of a systematic and general construction, unlike 
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FIG. 4: (a) A projection of the 3D cubic lattice onto the 
2D plane, (b) An example of a curve C with a framing 
C' . The links along the string are labeled 1,2,3. The links 
which cross C' are the filled dots labeled by 3, 4, 5, 6, 7. The 
corresponding twisted string operator W (C) is given by 
L+L 2 "L+(-l) 



VL Z ,+Li + Li+Li 



previous examples of emergent photons and fcrmions [7, 
12, 13], which were in some sense discovered by accident. 



A. The twisted string operator 

In order to motivate the twisted rotor model, we first 
discuss the "twisted string operator." 

We begin with the original "untwisted" rotor model 
-ffrotor- Notice that the operator B p = L^L^ L^L^ is a 
special case of a general string operator W(C) that can 
be associated with any curve C — *i*2*3---*ti in the cubic 
lattice: 



W{C) 



L+LVL+. 

ll «2 *3 



(6) 



B p corresponds to the case where C is the boundary of 
a plaquette dp. 

This "untwisted" string operator has an important 
property. When C is a closed curve, W{C) commutes 
with all charge operators and all other closed string op- 
erators: 



r, W(C)] = [W(C'),W{C)} = 



(7) 



This property is essential for the mapping to lattice gauge 
theory (4) to hold. It implies that the W(C) can be si- 
multaneously diagonalized within a given charge sector. 
The simultaneous eigenstates can then be interpreted as 
different flux configurations, and the W(C) can be in- 
terpreted as Wilson loop operators measuring the flux 
through a curve C. 

It turns out that there is another string operator - 
which we call the "twisted string operator" - that also 
has this property. The twisted string operator W(C) is 
similar to the usual string operator W(C) except that 
it contains an additional phase factor that depends on 
rotors on the "legs" of C - that is, rotors on the links 
incident to vertices Ik along the curve C. 

To define W(C), one must first choose a projection 
of the 3D lattice onto a 2D plane. Any projection will 
work; in this paper, we choose the particular projection 
illustrated in Fig. 4(a). The twisted string operator is 
then defined by 
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FIG. 5: A picture of the B p term in the twisted rotor model 
(9). Just like any other twisted string operator, the term 
B p = L+LzL+L±(-l) L i +L *+ L i+ L z+ L i+ L i+ L «+ L ™ acts on 
the four rotors, labeled by 1, 2, 3, 4, along the boundary of the 
plaquette p, and the 8 rotors, labeled by 1,4,5,6,7,8,9,10, 
that cross the framing curve. 



Here, n\ is the number of times that the link i crosses the 
curve C, where C is some "framing" of C. A "framing" 
is a curve C drawn next to C, but not exactly on it (see 
Fig. 4(b)). The choice of framing, C, is not unique, so 
whenever we discuss W, we will need to specify a partic- 
ular choice of framing. 

One can check directly that both of the equalities in 
(7) hold for the twisted string W . Alternatively, one can 
refer to the general argument given in Ref. [5]. 



The twisted rotor model 



The twisted model is obtained by modifying the B p 
term in (1). Instead of using a B p term based on the 



usual string operator, that is, 



Bp — 



W(dp), 



we use a 



W(C) = (LlL7. ..)(-!) 



1L: 



(8) 



modified B p term based on the twisted string operator: 
Bp = W(dp). We use a framing obtained by taking the 
plaquette boundary dp and shifting it up and to the left 
(see Fig. 5). The result is the twisted rotor Hamiltonian 

H% tm = VY,Q] + J^ L i) 2 - £s(Bp + h.c.) (9) 

lip 

logs of I 

Bp = L+L^L+Ll (-l) L l+ L l+ L l+ L i+ L 7+ L l+ L 9+Ll 

where the above explicit definition of B p applies to pla- 
quettes in the xz plane (see Fig. 5). (The definition for 
plaquettes in the xy and yz planes is similar and can 
be obtained using B p = W(dp) together with the above 
framing convention). 

Note that that (9) is very similar to the original rotor 
model (1) (and hence also very similar to lattice gauge 
theory) . The B p term alternately increases and decreases 
L z along the plaquette boundary dp, just like the B p 
term in the model (1). The only difference is that here the 
amplitude for this process is not always 1. The amplitude 
is ±1, depending on the rotors on the legs of the plaquette 
boundary. 
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C. Physical properties of the twisted rotor model 

The twisted rotor model can be analyzed in the same 
way as the untwisted model (Sec. II B). The charge oper- 
ators Qi commute with the Hamiltonian, so all the states 
are associated with a charge configuration {qi}. Low en- 
ergy states have charge-0 and can be thought of as string- 
net states. They can be enumerated in the same way as in 
the untwisted model. Just as the untwisted closed-string 
states could be obtained by repeatedly applying different 
Bp's to the vacuum state (the state with Lf — 0), the 
twisted closed-string states can be obtained by repeat- 
edly applying different Bp's to the vacuum state. Just 
as before, the J term can be interpreted as a string ten- 
sion, and the g term as a string kinetic energy. When 
g ^> J, string-net condensation occurs and the resulting 
ground state has two types of excitations: low energy 
string collective motions, and higher energy defects in 
the condensate. 

More quantitatively, one can show that the low en- 
ergy physics of the twisted rotor model is described by 
a compact U(l) gauge field coupled to infinitely mas- 
sive charges - just like the untwisted model. In fact, the 
twisted and untwisted rotor Hamiltonians (9), (1) are 
completely equivalent and can be mapped onto one an- 
other. 

To see this, note that the operators Qi and B p com- 
mute with each other and can therefore be simultane- 
ously diagonalized. Let the simultaneous eigenstates be 
denoted \{qi, p }, tw) where 

Qi\{qi,0 p },tw) = qi\{qi,0 p }, tw) 
B p \{qi, p }, tw) = e a '\{qi,e p },tw) 

and the qi, p are integers and real numbers, respectively. 
The states \{qi, P }, tw) form a complete basis for the 
twisted rotor model. Similarly, we can simultaneously 
diagonalize Qi and B p , to construct a complete basis 
\{li,0p}) f° r the untwisted model. 

Note that the matrix elements of B p , Qj in the twisted 
basis are the same as those of B p , Qi in the untwisted ba- 
sis. Furthermore the two sets of operators, (Qi, B p , Lf) 
and (Qi, Bp, L^), satisfy the same commutation rela- 
tions. Thus, Lj must also have the same matrix elements 
in the two models. We conclude that the two Hamiltoni- 
ans, (9) and (1), have exactly the same matrix elements 
(in their respective bases). The two models are therefore 
equivalent. 



D. Adding charge dynamics 

As before, we would like to modify the Hamiltonian 
-Orator s0 that the point defect excitations (or "charges" ) 
have a finite mass. To accomplish this, we need to add a 
term -ff^op tna ^ doesn't commute with Qj. Such a term 
will make the charges hop from site to site. 

We expect that any hopping term will give rise to the 
same qualitative behavior (e.g. same particle statistics). 
However, the resulting Hamiltonian is easier to analyze if 




FIG. 6: A picture of the hopping term HjZp = -tYU&t + 
Lt). The Lf are defined as twisted string operators W(i), 
with the framing curve drawn just below and to the right 
of i, as shown above. The result is Lf — L3 ( — l) L s+ i 6 
when i is the link labeled 3 pointing in the x direction, 
Lf = L* (~l) i?+L 8, when i is the link labeled 2 pointing 
in the y direction, and Lf — Lf( — l) L3+L4 when i is the link 
labeled 1 pointing in the z direction. 



we pick a hopping term which commutes with B p . While 
in the untwisted case the simplest hopping term hap- 
pened to have this property, in the case at hand we have 
to work a little harder. 

We notice that the untwisted hopping term ifhop is 
simply — t(W(i)+W(iy) where W(i) is the shortest open 
string operator (the string consists of a single link i). 
By analogy we choose the twisted hopping term to be 
— t(W(i) + W(iy) with a framing convention given by 
drawing the framing curve just down and to the right of 
the link i. The result is (see Fig. 6): 

i 

where 

(L±(-1) £ !+ £ S, if «|| a; 
Lt = { Lf(-l) L 'r+ L '*, if i || y (11) 
U±(-1) L 3+ L S if i|| x 

(One can check that Lf does in fact commute with B p 
- either by explicit calculation or using the general prop- 
erties of twisted string operators). [5] 

Note that the term flJJL causes the charges to hop 
from sites to neighboring sites. The low energy physics of 
(i?£rt 0r + #ho P ) i s therefore described by a compact [7(1 ) 
gauge field coupled to charges with finite mass. However, 
these charges are not Higgs bosons as in the untwisted 
case. Indeed, with the addition of charge dynamics, the 
two Hamiltonians (H^ tol . + H^ p ) , (H mtol +H hop ), are no 
longer equivalent. The two models are only equivalent in 
the special case where t = and the charges have infinite 
mass. As soon as the charges can hop, the mapping be- 
tween the two models breaks down, as does the mapping 
to lattice gauge theory. The reason for this is that the 
charged particles in (-ff*otor + -^hop) are n °t bosons, but 
fermions - as we will see in the next section. 
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E. The statistics of the U(l) charges 

In this section, we compute the statistics of the charges 
in the twisted rotor model. We use the same technique 
as in Sec. HE. 

As before, we perform the computation in the case 
J = 0. Since commutes with the Hamiltonian 

(-ff'^or + ffhop)> we can divide the Hilbert space into 
different sectors {9 P } corresponding to different flux con- 
figurations £>p W = e l6p . The states within each sector 
can be labeled by their charge configuration {qi}. The 
action of the Hamiltonian (H^ toI + H^ p ) on a charge 
configuration {qj} is simple. The first part of the Hamil- 
tonian, i?*otor> doesn't affect the charge configuration at 
all, while the second part, H^ p acts m tw0 ways: it ei- 
ther creates two charges at neighboring sites, or it makes 
a charge hop from one site to another. 

Thus, within each sector, (-ff'otor + ^ho P ) * s simply a 
hopping Hamiltonian on the cubic lattice. The hopping 
operators Lf make the charges hop from site I to site J, 
where 7, J are the two endpoints of i. To compute the 
statistics of the charges, we use the statistical hopping 
operator algebra [11], just as in the untwisted case. We 
find that the hopping operators satisfy the same algebra 
but with a — sign. That is, 

LtLjL+ = (-l)-L+LjL+ (12) 

for any i, j, k incident to some vertex J. This sign makes 
all the difference. The above algebra is not the bosonic 
hopping operator algebra, but rather the fermionic hop- 
ping algebra. We conclude that in the twisted rotor 
model the charged particles are fermions. 

One might expect that, by analogy with the untwisted 
case, there is a mapping to fermionic lattice gauge theory. 
However, no such mapping exists. The reason is that, at 
the lattice scale, these fermions behave differently from 
the usual fermions in a tight binding model. 

There are two key differences. First, in a standard 
tight binding model of spinless fermions, each site can 
be occupied by or 1 fcrmion. However, in our case the 
charge Qi can be any integer from — oo to oo. Thus, a 
given site can be occupied by arbitrarily many fermions. 
Second, in a standard tight binding model, the fermions 
are totally noninteracting in the zero coupling limit J — 
0. However, in our case the the fermions do interact, even 
in this limit. For example, the hopping term i?hop allows 
positive and negatively charged fermions to annihilate 
each other if they occupy neighboring sites. 

Neither of these differences should affect the physics of 
the twisted rotor model in the low energy limit, if V is 
large. In that case, the low energy physics of the rotor 
model is equivalent to an insulator: it is described by two 
species of gapped fermions with opposite charges coupled 
to a U(l) gauge field. 

However, these microscopic differences could be impor- 
tant when V becomes small enough that a Fermi surface 
develops. Unfortunately, in order to construct massless 
Dirac fermions, we need to consider that limit. Therefore, 
in order to make our calculation well-controlled, we need 
to modify the twisted rotor model so that the fermions 



are equivalent to tight binding model fermions, even at 
the lattice scale. This is the subject of the next section. 



IV. ENGINEERING MASSLESS DIRAC 
FERMIONS 

In this section, we modify the twisted rotor model 
(H^ toI + H^ p ) so that the fermions are exactly equiva- 
lent to tight binding model fermions. We then show that 
this more well-controlled model can give rise to mass- 
less Dirac fermions. It can also give rise to Fermi liquid 
behavior. 



A. The model 

There are many ways to make the fermions in (H^ tOT + 
flhop) behave like the usual fermions from a tight binding 
model. Here we will only describe one of them. First, 
we divide the vertices of the cubic lattice into 3 classes: 
"eee vertices", "ooo vertices" and "mixed vertices." A 
vertex is called an eee vertex if the three components 
I = (I Xl I y ,I z ) are all even, an ooo vertex if the three 
components J = {I x ,I y ,I z ) are all odd, and a mixed 
vertex otherwise. We then modify the V term in (1) to 

\ eee ooo mixed / 

Throughout our discussion we will take the V — > +oo 
limit. In this limit, Qi must be 0, 1 on the eee vertices, 
0, —1 on the ooo vertices, and on the mixed vertices. In 
other words, the U(l) charges can only live on the eee or 
ooo vertices, and these two sublattices can only contain 
unit positive and negative charges respectively. 

Next, we modify the hopping Hamiltonian to 

-t ]T (LiLT+LTL+) (13) 

(ij),coll. 

where the sum runs over pairs of neighboring links (ij) 
which arc collinear and can be joined together to form a 
single line segment of length 2. This hopping term has 
the property that charges only move within the eee and 
ooo sublattices. 

Putting this all together, and including a term which 
will correspond to a chemical potential /x, we arrive at 
the following rotor model: 

#t.b - J^K) 2 -g^2(B P + h.c.) - 2^ L l 

i p i 

+ v (E Q^Q* - 1) + E QilQ* + + E & ) 

\ eee ooo mixed / 

"* E (LtLj+L7L+)) (14) 

(ii), coll. 
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B. Solving the model 

We begin by considering the case where t, J = 0. In 
this case -fft.b is exactly soluble. Indeed, the basis states 
$p}> tw) are eigenstates of H t ^ with an energy 

£ = W^ w ( gf -l) + ^ gi ( 9l + l) + Yl 1i\ 

\ eee ooo mixed / 

-m£(-1)V-2 5 £;cos(0p) (15) 
i P 

As we mentioned earlier, we will take the U — > oo limit. 
In this limit, the only low energy states are those with 
qi = 0, 1 in the eee sublattice, qi = 0, —1 in the ooo 
sublattice and qi — everywhere else. If we restrict to 
this low energy subspace, the energy E is simply given 

by 

E = -fj,N - 2g cos(e p ) (16) 
p 

where TV = J^eee I 1 ~ Tlooo I 1 i s tne total number of 
charged particles. 

Now, consider the case t ^ 0. The t term generates 
nearest neighbor hopping in the eee and ooo sublattices. 
If we let P denote the projection onto the low energy 
subspace, then the action of the hopping Hamiltonian 
(13) within this subspace can be written as 

Y (*T(I+2»/) + h.C.) + X! (tl(I+2") + h - C -) 
eee.v ooo,v 

where u runs over the unit vectors x, y, z and the tj(j +2u ) 
are hopping operators that make the particles hop from 
site I + 2is to site I. They are defined by 

ti(i+2„) = -t ■ PLflJP (17) 

where i is the link connecting I, I + u, j is the link 
connecting I + u, I + 2v, and we use the upper (lower) 
signs for the eee (ooo) sublattices. 

The low energy effective Hamiltonian can then be writ- 
ten as 

H t .b, eff = ^ {tl(I+2u) + h.C.) 

eee,f 

+ (*-T(J+2f) + h - C -) 

ooo.v 

- fiN- 2gJ2 cos ( d P ) ( 18 ) 
v 

This Hamiltonian describes two types of hardcore 
charged particles with charges ±1 hopping on two dif- 
ferent sublattices. We can completely characterize this 
Hamiltonian by investigating two physical properties: 
the statistics of the charges, and the gauge flux through 
each plaquette p. 

To determine the statistics of the charges, we use the 
statistical hopping operator algebra. [11] We notice that 
the hopping operators satisfy the relation 

tlLtKltlJ = (-1) • hjtKltlL (19) 



for any J,K,L adjacent to J. According to [11], this is 
the fermionic hopping operator algebra, so the particles 
arc fcrmions. 

To determine the flux that the particles see, we com- 
pute the product of hopping operators around a 2 x 2 
plaquette (IJKL) in the eee or ooo sublattice. We find 
that (up to some signs having to do with orientation con- 
ventions), 

tljtjKtKLtLI (20) 

= (-t) 4 ■ e^^uKL) i0 v ni (i _ nj )(i _ n K )(l - n L ) 

Here the sum on the right hand side runs over the 4 
plaquettes in the original lattice that are contained in 
the (doubled) plaquette {IJKL). Also, nj = (— l) I Qi 
denotes the occupation number of site I. This relation 
implies that the particles see a flux P through each pla- 
quette p. 

The hopping operators tjj are completely character- 
ized by the two algebraic relations (19), (20). Any col- 
lection of hardcore hopping operators that satisfy these 
relations will give rise to a Hamiltonian equivalent to 
(18). A particularly simple collection of hopping oper- 
ators satisfying these relations can be constructed from 
fermionic operators cj with 

{c/,cj} = 0, {c I ,c] T } = 5 IJ 

The hopping operators are given by 

ii(i+2u) = ~t ■ e ±li > A ^+»<> +A v+^*+^c\c I+2 „ (21) 

where the upper (lower) sign applies to the eee (ooo) 
sublattice and e lAlJ are phases defined on each link sat- 
isfying 

e iA 'iJ» e iA '.'«-i e tA ' 3 i2e lAl 2'i =e T. ve c l9 p 

for any closed curve C. In other words, A is a vector 
potential for the flux configuration 8 p . 

Rewriting the low energy effective Hamiltonian in 
terms of these fermionic operators gives 

Ht.b, off' 

= Y ( _ * " e i ( Al <. I +^ +A V+»W+ 2 ^ ) c\c I+2 u + h.c.) 

eee,v 

+ H ( _ * ' e-* (Aj ( J +">- M < I +">< I + 2 "> ) c\c I+2v + h.c.) 

ooo,is 

~M X! C J C/ ~ ^ ^2 °i Cl ~ 2g ^2 C0S (8 P ) ( 22 ) 

eee ooo p 

Note that this is nothing other than the standard tight 
binding Hamiltonian for fermions hopping in fixed flux 
configurations 8 p . If wc now allow J ^ 0, then the flux 
configurations acquire dynamics. The Hamiltonian (22) 
becomes modified by the addition of an electric energy 
term JJ2{u) Ejj, where Eu is canonically conjugate 
to Aij: \Ajj , Eu] = i. This term is precisely what 
one needs for compact lattice U(l) gauge theory. We 
conclude that in the general case J, t ^ 0, the twisted 
rotor model (14) is mathematically equivalent to a tight 
binding model with ±1 charged fermions coupled to a 
compact U(l) gauge field. 
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C. Fermi liquid behavior 

The Hamiltonian i7t.b,cff naturally gives rise to Fermi 
liquid behavior. This is easiest to understand in the 
regime g 3> |J|,|i|. In this case, the minimum energy 
flux configuration is the configuration 6 P = 0. The flux 
6 P fluctuates about this minima giving rise to photon- like 
excitations. The fcrmion dispersion in the background 
flux configuration is given by 

E(k) = -2t(cos(2k x a) + cos(2k y a) + cos(2k z a))-fi (23) 

The ground state is obtained by filling all the negative 
energy levels. Clearly, the model behaves differently for 
different values of chemical potential fi. If /i < —6t, 
then there are no negative energy levels. All the energy 
levels are empty in the ground state, and the system is 
an insulator. In this case, the rotor model (14) describes 
gapped fermions with a gap A = — 6t — \i coupled to 
a compact 17(1) gauge field. The fermions come in two 
species (corresponding to the two sublattices, eee and 
ooo) with charges +1 and —1, respectively. Hence, just 
like the original twisted rotor model (H^ toI + ), the 
physics is equivalent to that of a 17(1) gauge field coupled 
to two species of gapped fermions with opposite charges. 

On the other hand, if fi > — 6t, then the result is a 
Fermi liquid - or more precisely two Fermi liquids corre- 
sponding to the two species of fermions. These two Fermi 
liquids arc both coupled to a compact 17(1) gauge field. 

D. Massless Dirac fermions 

The model (22) can also give rise to massless Dirac 
fermions. Indeed, massless Dirac fermions can occur 
whenever the fermion band structure has nodal points. 
This happens naturally in a 7r-flux configuration where 
each plaquette contains a flux of it. 

In our case, the relevant plaquettes are not plaquettes 
p in the original cubic lattice, but rather "doubled" 2x2 
plaquettes in the eee and ooo sublattices. One way to 
ensure that the flux through these "doubled" plaquettes 
is 7r, is to introduce a spatial dependence into the param- 
eter t. That is replace t by t ■ \u where 

Xi,i+2x =1) XiJ+2 V =(-l) ix/ ' 2 , 

Xi.i+2* =(-l)( 7 -+ J »>/ 2 (24) 

on the eee sublattice and 

XU+2X =1, Xi,i+2y =(-l)7*-D/ 2 , 

Xi,i+2z =(-l)( J -- 1 + J »- 1 )/ 2 (25) 

on the ooo sublattice. The result is that (20) becomes 
modified so that each fermion sees an effective flux 
°f Tlpe(iJKL) ®p + n through each doubled plaquette 
(IJKL). 

The minimum energy flux configuration is still 9 P = 0, 
but now the fermions see an effective flux of ir through 
each doubled plaquette. The fermion dispersion in this 



background is given by 2 pairs of degenerate bands with 
energies 



E(k) — ±ty cos 2 (2fca,a) + cos 2 (2k y a) + cos 2 (2fc z a) - /j, 

(26) 

and Brilluoin zone, (— 7r/4a, 7r/4a) x (— 7r/4a, ■n/Aa) x 
(— n/2a, 7r/2a). The ground state of H t .b, off is again ob- 
tained by filling all the levels with negative energies. If 
/i = 0, then the Fermi surface consists of 2 nodal points 
with linear dispersion - at k = (ir/Aa, it /4a, ±7r/4a). The 
low energy fermionic excitations are concentrated near 
these 2 nodes. Together, the two nodes give rise to one 
four-component Dirac fcrmion in the continuum limit. 
Since the bands are doubly degenerate, and there are 
two species of fermions (corresponding to the eee and 
ooo sublattice), there are a total of 4 massless four com- 
ponent Dirac fermions. 

We conclude that the low energy physics of the rotor 
model H t .b with n = 0, V — ► oo, g ~> J,t, and t given 
by the above formula, is equivalent to QED with four 
species of massless Dirac fermions. Thus, 7/t.b is a lo- 
cal bosonic model that gives rise to both photons and 
massless fermions. It can be viewed as a realization of a 
quantum ether - an ether that produces not only light, 
but also electrons! 

We would like to point out that the speed of light c ~ 
\fg~J a and the speed of the massless Dirac fermion c e ~ ta 
are not the same in our rotor model. However, we can 
tune the value of t to make c = c e . In this case the low 
energy effective theory will have Lorentz invariance. The 
reason that we do not get a Lorentz invariant low energy 
effective theory naturally is that our approach is based on 
Hamiltonian formalism where space and time are treated 
very differently. If one uses a path integral formalism 
and discretizes space and time in the same way, Lorentz 
invariance can emerge naturally in a low energy effective 
theory. 

Another difficulty with the above model is that the 
emergent electrons are massless. However, this prob- 
lem can also be overcome: when bosonic models contain 
emergent non-Abelian gauge bosons, the models can pro- 
duce electrons with a finite mass and such a mass can be 
much less than the cut-off scale without fine tuning any 
parameters. [13] 



V. CONCLUSION 

In this paper we have explicitly demonstrated the 
string-net picture of (3+l)D emerging gauge bosons and 
fermions. We presented a rotor model on the cubic lat- 
tice with a string-net condensed ground state, and ex- 
citations which behave just like massless photons and 
massless Dirac fermions. We saw that the model was 
closely related to 17(1) lattice gauge theory coupled to a 
Higgs field. The only real difference was a phase factor 
or "twist" in the Hamiltonian. Surprisingly, this "twist" 
was enough to turn the Higgs boson into a fermion. 

This model is part of a much more general construc- 
tion [5] that can produce gauge bosons with any gauge 



10 



group. In particular, one should be able to construct a 
model that gives rise to SU(3) gauge bosons and massless 
Dirac fermions - that is, QCD. These models could have 
applications to lattice gauge theory. They may make it 
possible to simulate QED with electrons and QCD with 
quarks without ever introducing quark or electron de- 
grees of freedom on the sites. Only the link variables 
(representing the gauge field) would be necessary in such 
a simulation. 

From a high energy point of view, the string-net pic- 
ture of the vacuum is quite appealing. It explains why 
the standard model looks the way it does - that is, why 
nature chooses such peculiar things as gauge bosons and 
fermions to describe itself. In addition, it unifies the mys- 
terious gauge symmetries and anticommuting fields into 
a single underlying structure: a string-net condensate. 

But can we actually construct a string-condensed local 
bosonic model that produces the entire standard model? 
We are close, but not quite there. In terms of elementary 
particles we can produce photons, gluons, leptons and 
quarks, but we do not know how to produce neutrinos or 
SU(2) gauge bosons. 



The problem with the neutrinos and the SU (2) gauge 
bosons is the famous chiral-fermion problem. [19] Neu- 
trinos are chiral fermions and the SU(2) gauge bosons 
couple chirally to other fermions. At the moment, we do 
not know how to obtain chiral fermions and chiral gauge 
theories from any local lattice model, much less a local 
bosonic model. 

From a condensed matter point of view, the above 
model shows how a simple spin system can give rise to 
emergent fermions and gauge bosons. Extended objects 
can give rise to both of these phenomena, as long as their 
condensate has an appropriate "twist." This may pro- 
vide intuition in the search for new and exotic phases of 
matter - phases beyond the scope of Landay's theory of 
symmetry breaking. The discovery of a real material with 
string-net condensation would represent a breakthrough, 
particularly a material containing excitations which be- 
have just like the photons and electrons in our vacuum. 

This research is supported by NSF Grant No. DMR- 
04-33632 and by NSF-MRSEC Grant No. DMR-02- 
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